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I.  INTRODUCTION: 


The  nitromethyl  radical  has  been  postulated  to  be  involved  in  a  number 
of  different  phenomena*  In  a  study*  of  nitromethylation  of  aromatic 
compounds  it  was  determined  that  the  radical  had  appreciable  electrophilic 
character*  The  only  products  formed  were  through  reactions  at  the  carbon 
suggesting  that  it  was  the  more  reactive  center.  Interestingly  in  the 
reaction  with  manganese  (III)  acetate  and  nitrome thane  the  authors* 
suggested  that  it  might  be  the  aci  form  (CH^ONO)  that  underwent  oxidation 
followed  by  tautomerization* 

The  radical  could  be  produced  in  the  following  reaction: 

ch3  +  ch3no2  +  ch4  +  ch2no2 

which  has  been  suggested2  to  explain  the  appearance  of  methane  as  a  major 
product  in  the  decompostion  of  nitromethane.  Once  formed  the  nitromethyl 

radical  is  believed  to  effect  the  thermal  decomposition  of  nitromethane  in 

3  .  ... 

condensed  phase  systems  .  Such  systems  include  propellant  ignition  and 

explosive  initiation*  Nitromethyl  radical  may  even  play  a  role  in  the 

chemistry  of  the  upper  troposphere  through  the  degradation  of  nitromethane 

.  «  4 

formed  from  the  possible  reaction  of  CH3  and  NO^. 

Very  recently  the  vibrational  spectrum  of  CH^NO^  was  reported  by 

4  z  1 

Jacox  who  codeposited  F  atoms  and  nitromethane  zn  an  Argon  matrix  at  14 

K  and  analyzed  the  changes  in  the  IR  spectra*  Different  isotopes  of 

nitromethane  were  used  in  the  study  to  help  assign  the  bands* 

The  heat  of  formation  of  gaseous  nitromethyl  radical  has  been 

estiraatm  to  be  37*3  kcal/mol* 

The  ESR  spectrum  of  the  nitromethyl  radical  has  been  interpreted^  as 
supporting  a  nonplanar  geometry.  However,  the  only  reported  theoretical 
calculation  on  CH2N02  was  an  INDO  calculation^  which  predicted  a 
planar  structure. 

The  motivation  for  this  study  was  not  only  to  gain  an  understanding  of 
this  particular  system  as  a  reactive  intermediate,  but  also  to  view  the 
system  as  a  model  for  the  radicals  produced  from  high  energy  species  such 

g 

as  TNT  (1,3,5-trinitrotoulene)  where  radical  intermediates  as  1  or  2 
might  occur. 


The  two  valence  structures  which  can  be  written  for  planar 


CH2N02  are: 


1 


3  4 


A  theoretical  study  on  radicals  of  this  nature  involving 
structures  with  comparable  energies  but  inequivalent  assignments  of 
ir  bonds  and  the  lone  electron  suggest  that  accurate  results  require 
substantial  effort  to  correctly  predict  their  relative  energy*  Indeed 
this  could  be  predicted  by  the  sophisticated  basis  sets  needed*®  to 
predict  the  states  of  NO^. 

II.  OBJECTIVES  OF  THE  RESEARCH  EFFORT: 


The  object  of  this  work  was  to  study  the  types  of  radicals  that 
could  be  generated  by  thermal  or  photolytic  decomposition  of  high 
energy  materials.  Specifically  the  radicals  produced  from 
polynitrated  systems  as  TNT  or  HMX.  To  do  this  one  must  determine 
which  is  the  most  economical  method  that  will  produce  results  of 
sufficient  quality  or  accuracy  to  answer  the  particular  chemical 
problem*  If  empirical  methods  are  not  sufficiently  accurate  and  more 
sophisticated  methods  too  costly,  one  may  ask  if  there  are  trends 
which  could  still  allow  the  empirical  methods  to  be  used  predictively 
with  some  confidence.  The  objectives  of  this  study  were  as  follows: 

(1)  To  study  the  electronic  states  of  the  CH^NO^  radical  with 
the  semiempirical  method  MNDO  to  determine  the  ground  state  and 
electronic  configuration. 

(2)  To  study  the  problem  using  ab  initio  methods  with  geometry 
optimization  and  correlation. 

(3)  To  determine  the  extent  that  other  configurations  play  in 
determining  the  properties  of  the  electronic  states  by  using  a  MCSCF 
method. 
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III.  RESULTS  AND  DISCUSSION: 


The  first  method  used  to  study  the  electronic  states  of 
was  the  MNDO  method,  a  seraiempirical  method  based  on  the  NDDO 
approximation**  and  very  successfully  parameterized  by  Dewar  and 
Thiel  .  All  calculations  were  carried  out  using  the  unrestricted 
SCF  formalism. 

The  heats  of  formation  and  geometries  for  CL  species  are  given 

2 

in  Table  1.  The  most  stable  state  is  the  planar  B,  followed  very 

2  . 

closely  by  the  staggered  B2 .  These  two  states  simply  correspond 

to  rotation  of  the  singly  occupied  p  orbital  on  carbon.  The  other 

states  were  obtained  by  starting  with  a  bond  order  generated  from  an 

orbital  guess  in  which  selected  orbitals  were  occupied.  During  the 

iteration  to  self  consistency  the  switched  orbital  was  monitored  to 

insure  that  it  remained  occupied.  Unfortunately  in  many  instances  the 

wavefunction  lost  the  initial  C2v  symmetry  and  it  proved  to  be 

impossible  to  obtain  those  states. 

2 

The  planar  A9  solution  was  found  to  be  18.3  kcal/mol  higher 

2  2 

in  energy  than  the  planar  B1  state.  Roughly,  these  two  states 

correspond  to  valence  structures  3  and  4.  The  unrestricted  solution 

2 

obtained  for  A^  was  characterized  by  three  singly  occupied 
orbitals;  a  b^  and  a2  orbital  of  alpha  spin  and  a  b^  orbital  of 
beta  spin. 

The  assumption  of  symmetry  for  radical  seems 

intuitively  obvious  but  due  to  the  interaction  of  closely  related 
electronic  state  this  assumption  may  be  wrong.  It  has  previously  been 
noted*^a,k  that  under  certain  conditions  electronic  states  in 
symmetry  may  mix  by  an  asymmetric  distortion  and  lead  to  a  lower 
energy  in  symmetry.  The  conditions  which  must  be  met  are: 

1)  the  two  states  in  C ^  symmetry  must  transform  as  the 

same  irreducible  representat ion  in  the  C  point  group; 

s 

2)  there  must  be  a  surface  crossing  at  some  nearby  C2v 

geometry; 

3)  the  leading  configurations  must  differ  by  only  one  spin 

orbital. 

The  ^A-  and  states  of  NO-  meet  the  above  conditions 

1  4  A  1  »  “ 

and  it  is  found1  a’  that  an  asymmetric  distortion  does  lead  to  a 


3 


lower  energy. 

If  we  look  at  the  C^v  point  group  table  we  see  that  the  C2v 
symmetry  can  be  reduced  to  Cg  by  the  loss  of  the  sigma  (yz)  plane  or 
the  sigma  (xz)  plane  (J>) 


If  the  xz  symmetry  plane  is  lost  then  the  A2  and  states 

transform  as  A"  while  A^  and  B2  states  transform  as  A1  •  If  the  yz 
symmetry  plane  is  lost  then  and  B2  transform  as  A"  while  the 

A^  and  B^  transform  as  A*. 

The  MNDO  wavefunct ions  were  tested  for  stability  under  Cg 
distortions  and  the  results  are  presented  in  Table  2  and  3. 

Table  1  MNDO  Heats  of  Formation  (kcal/mol)  and  Geometries 
of  radicals  in  C2v  symmetry 


State 

AH, 

CN 

NO 

CH 

ONH 

HCN 

li 

i 

35.2 

1.476 

1.214 

1.084 

119.3 

118.3 

4 

53.5 

1.414 

1.266 

1.084 

122.3 

119.0 

B2 

35.6 

1.482 

1.219 

1.082 

119.4 

117.6 

2Ac,d 

A1 

63.1 

1.454 

1.264 

1.084 

121.5 

118.2 

a)  Distances  are  given  in  Angstroms  and  angles  are  given  in  degrees. 

b)  Planar  geometry. 

c)  Staggered  geometry. 

2 

d)  The  A^  state  has  b^ ,  ,  and  b^  orbitals  singly  occupied. 

Table  2  MNDO  Heats  of  Formation  (kcal/mol)  on  distortion  from 
C2v  Symmetry  to  Cg  Symmetry 

Element  of  Symmetry  Symmetry 

»  a  b 

symmetry  kept  in  C2y  in  Cg  AH^  AH^ 

sigma(yz)  2A,,  2B.  2A"  33. 8C  (1.4)  71.0 


sigma(xz) 


33. 8C  (1.8) 


d 


2a  2r 

2.  2„ 

Al*  B1 


61.6  (1.5) 


a)  Planar  geometry*  Number  in  parenthesis  refers  to  lowering  in 
kcal/mol  with  respect  to  the  lower  of  the  two  C2v  states. 

b)  Staggered  geometry.  Number  in  parenthesis  refers  to  lowering  in 
kcal/mol  with  respect  to  the  lower  of  the  two  C^v  states. 

c)  The  wavefunction  exhibited  only  symmetry. 

d)  A  distortion  from  C2V  symmetry  was  not  energetically  favorable. 


Table  3  MNDO  Geometries 
of  CH^NO^  radicals  in  Cg  symmetry 


Geometry  AHf  CN  NO  CH  ONC  HCN  Dihedral 

6  33.8  1.466  1.221  1.083  120.5  118.3 

1.227  1.083  119.1  118.5 

6a  61.6  1.465  1.268  1.084  118.0  120.2  72.8 

1.082  116.4 

7  71.0  1.452  1.226  1.084  125.3  118.0  88.9 

1.306  116.6 

a)  When  two  values  for  a  variable  are  given,  the  first  value 
corresponds  to  the  bond  distance  or  bond  angle  on  the  side  of  the  CN 
bond  denoted  'a*  in  figures  6^,  6a ,  and  7^  while  the  second  value 
corresponds  to  the  side  denoted  *  b 1 . 


The  global  minimum  was  a  distorted  planar  structure  6  but  since 

the  wavefunction  did  not  possess  C  symmetry,  the  state  could  not  be 

s  2  2 

definitively  assigned.  Although  the  staggered  A^  and 

solutions  could  not  be  found,  the  MNDO  value  of  71.0  kcal/mol  for  the 

heat  of  formation  of  the  distorted  structure  6a  is  a  lower  bound  for 

the  lower  of  the  states.  When  the  sigma  (yz)  plane  of  symmetry 

2  2  2 
is  lost  the  A.  and  B,  states  interact  and  lead  to  a  A1 

11  2 
state  2  1«5  kcal/mol  lower  in  energy  than  the  A^  state. 


6  6a  7 
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Since  the  relative  energies  of  NO^  and  related  compounds  are 

particularly  sensitive  to  the  basis  set^,  ab  initio  calculations 

13 

were  performed  on  the  CH^NO^  system  •  Geometries  were 

completely  optimized  within  either  the  C~  or  C  point  group  using 
14  s 

the  3-21G  basis  and  single  point  calculations  were  then  made 
including  correlation  at  the  UMP3  level^^.  The  energies  relative  to 


2 

the  staggered  are 

presented  in  Table  4S 

absolute  energies  in 

Table  3,  and  geometries 

in  Table  6 

•  From  Table  4  one  can  determine 

that  the 

2  2 
%  "  Bi 

(planar) 

separations  calculated  at  the 

UHF/3-21G, 

UHF/6-31G,  UMP2/6-31G, 

and  UMP3/6- 

-31G  are  respectively; 

-17.7,  -4.6 

,  49.4,  and  31 

.0  kcal/mol, 

This  compares  to  the 

MNDO  value 

of  18.3  kcal/mol  from  Table  1. 

Table  < 

4  Relative  Energies  (kcal/mol)  for  Different  States  of 

CH^NO^  Based  on  Single  Configuration 

Calculation 

UMP3/6-31G 

UMP2/6-31G 

UHF/6-31G 

UHF/3-21G 

o 

Planar  A~ 

24.3 

45.1 

-12.7 

-26.5 

\ 

-6.7 

-4.3 

-8.1 

-8.8 

\ 

16.5 

29.8 

9.2 

0.3 

\ 

122.9 

2a" 

5.7 

23.3 

-22.8 

-32.6 

2a* 

37.5 

57.9 

9.0 

-16.0 

Stagg.  Zk 

75.4 

87.7 

58.2 

45.4 

\ 

110.4 

133.0 

76.0 

53.1 

\ 

58.1 

71.6 

38.6 

30.9 

\ 

0 

0 

0 

0 

\ 

80.7 

2A  " 

52.5 

73.7 

15.5 

0.9 

2a  ' 

18.0 

36.9 

-14.4 

-24.6 

Table  5  Total  Energies  (Hartrees)  for  Different  States  of 
C**2N02  ®ased  on  Single  Configuration  Calculation 


State 

Planar 

:  config.3 

UMP3/6-31G 

UMP2/6-31G 

UHF/6-31G 

UHF/3-21G 

16/3/2/10 

243.33871 

243.35832 

242.89561 

241.62565 

»2 

16/4/2/9 

243.30166 

243.30401 

242.86808 

241.61118 

v 

£ 

16/4/1/10 

16/5/1/9 

243.28929 

243.27959 

242.90298 

241.65386 

241.41585 

2A"b 

26/5 

243.31891 

243.31428 

242.91899 

241.66365 

(  6.1)c 

2A,d  25/6 

Staggered 

243.26825 

243.25916 

242.86835 

241.63718 

(16.3)c 

16/5/2/8 

243.23547 

243.23744 

242.82115 

241.56241 

% 

16/4/2/9 

243.32802 

243.35148 

242.88272 

241.61165 

\ 

15/5/2/9 

243.20787 

243.21167 

242.78989 

241.53924 

\ 

\ 

16/4/1/10 

16/5/1/9 

243.15204 

243.13955 

242.76159 

241.52701 

241.48298 

2A"C 

24/7 

243.24441 

243.23404 

242.85805 

241.61020 

(30.0)c 

2A'f 

25/6 

243.29926 

243.29263 

242.90567 

241.65085 

(24.6)c 

a)  The  notation  refers  to  the  number  of  electrons  in  the  a^,  b^, 

a«,  and  b^  orbitals  respectively.  The  notation  for  ^A1'  and 

A1  refers  to  the  number  of  electrons  in  a1  and  aM  orbitals 
respectively 

2  2 

b)  Involves  the  interaction  of  planar  states  and  A^  and 

the  loss  of  the  sigma  (xz)  plane  of  symmetry. 

c)  The  number  in  parenthesis  refers  to  the  lowering  in  kcal/mol  with 
respect  to  the  lower  of  the  two  interacting  C2v  states  for  the  3-2 1G 
energies. 

2  2 

d)  Involves  the  interaction  of  planar  states  B2  and  A^  and 

the  loss  of  the  sigma  (xz)  plane  of  symmetry. 

2  2 

e)  Involves  the  interaction  of  staggered  states  B^  and  A^ 

and  the  loss  of  the  sigma  (xz)  plane  of  symmetry. 

2  2 

f)  Involves  the  interaction  of  staggered  states  B^  and  A^ 

and  the  loss  of  the  sigma  (xz)  plane  of  symmetry. 
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Table  6  Ab  initio  UHF/3-21G  Geometries  of  CH^NO^  radicals 
in  Cjy  Symmetry  and  the  Spin  Squared  Values  (without  projection) 


2 

State  S 

CN 

NO 

CH 

ONC 

HCN 

Dihedral 

Planar 

*2 

1.59 

1.377 

1.360 

1.065 

122.5 

117.2 

0.79 

1.408 

1.257 

1.065 

117.1 

116.6 

2b2 

1.08 

1.281 

1.369 

1.066 

129.2 

118.3 

\ 

1.80 

1.384 

1.428 

1.064 

138.6 

116.1 

2A"a 

1.47 

1.391 

1.288 

1.065 

124.5 

117.1 

1.362 

1.064 

115.1 

117.1 

2A’2 * * 5 * * 8 

1.37 

1.341 

1.458 

1.065 

115.0 

118.9 

1.284 

1.066 

129.8 

116.5 

Stagg. 

A 

1.54 

1.336 

1.454 

1.068 

128.4 

117.2 

4 

1.84 

1.425 

1.370 

1.068 

115.2 

116.7 

2»i 

1.84 

1.413 

1.363 

1.070 

128.2 

117.8 

% 

0.76 

1.456 

1.241 

1.067 

117.0 

116.5 

2a. 

1.84 

1.402 

1.419 

1,071 

124.7 

118.0 

2a,,c 

1.76 

1.419 

1.472 

1.070 

111.8 

118.7 

90.6 

1.299 

128.9 

2a,c 

1.63 

1.431 

1.319 

1.070 

120.8 

118.0 

96*1 

1.328 

118.7 

a)  See 

Table 

3  and  Figure 

6. 

2 

b)  The  geometries  for  the  two  perpendicular  states  of  A£ 
symmetries  differ  in  the  orbital  occupancy*  The  first  entry  was 
obtained  with  3  open  shells,  and  the  second,  with  one  open  shell(see 
text) • 

c)  See  Table  3  and  Figure  8* 

2  2 

The  rotational  barrier  of  planar  to  staggered  is 

fairly  insensitive  to  the  level  of  calculation  (7-9  kcal/mol)*  The 

2  2 

rotational  barrier  of  planar  to  staggered  A^  is  also 

insensitive  to  correlation  if  the  same  configuration  is  used*  The 

barrier  at  the  UHF/6-31G  level  is  79*6  kcal/mol  compared  to  86*1 

kcal/mol  at  the  UMP3/6-31G  level.  The  high  barrier  reflects  the 

multiple  bond  character  of  the  ON  bond*  At  the 

single  determinant  level,  however,  the  [8at2b.5b9] la*  configuration 

5  2  *2  *2  ^  1  *1  1  16 

leads  to  a  higher  energy  than  the  papb^lajAb^lSapb *5b2  configuration 

The  energy  difference  between  the  two  solutions  varied  from  10-30 

kcal/mol  depending  on  level  of  calculation  (see  Table  4). 


2 

Of  course  Che  true  staggered  state  will  be  some  linear 

combination  of  these  two  states  plus  other  configurations  of  the 
correct  symmetry.  The  proper  treatment  of  such  states  must  be  based 
on  either  an  MCSCF  approach  or  a  full  configuration  interaction 
approach. 

2 

As  seen  in  Table  6  the  CN  bond  length  of  the  planar  B.  states 

2  1 

elongates  0.05  A  upon  twisting  to  the  staggered  B2  state.  This 
suggests  that  the  6.7  kc a 1/mol  required  for  rotation  (UMP3/6-31G)  may 
be  due  either  to  elimination  of  conjugation  in  the  planar  structure  or 
increased  repulsion  of  the  electron  in  the  singly  occupied  orbital 
with  the  oxygen  atoms  in  the  staggered  structure. 

It  must  be  remembered  that  this  is  not  an  electronic  transition 

between  states  but  rather  a  conformational  change.  The  states 
2  2 

B^  and  B^  are  antisymmetric  with  respect  to  C 2  rotation 

while  A2  state  is  symmetric. 

In  the  former  case  a  rotation  will  split  the  and 

2  2  ^  ^ 
planar  states  into  1  A”  and  2  A"  which  may  interact.  The  rotation 

continues  until  at  180  degrees  the  states  become  noninteracting  with 
,  .  2  2 

new  designations  of  B^  and  B^ .  The  interaction  could  lead 

to  a  weak  minimum  for  the  lower  state  or  a  slight  maximum  in  the 

higher  state.  However  in  this  paper  rotational  barrier  refers  to  the 

2 

energy  difference  in  the  two  C_  electronic  states.  For  A,. 

2  2^ 

rotation  the  same  applies  to  the  interaction  of  A~  and  A.  as 
2  2.  21 
1  A1  and  2  A1.  In  this  case  the  energy  separation  between  the  two 

states  make 8  mixing  unlikely. 

2 

The  shortest  CN  bond  occurs  in  the  planar  B.  state,  0.10  A 

2  ^ 

shorter  than  the  CN  bond  in  the  planar  A^  state  4  whose  valence 

structure  can  be  drawn  with  a  double  bond.  The  shorter  CN  bond  is  due 

to  removing  an  electron  (with  respect  to  the  A^  configuration) 

from  the  CN  antibonding  5b^  orbital  and  putting  it  in  the  CN  bonding 

la2  orbital.  The  loss  of  a  electron  from  the  5b2  orbital  also 

causes  the  CNO  bond  angle  to  increase  about  6  degrees  when  compared  to 
2 

the  planar  A^  geometry. 

The  spin  squared  values  for  the  3-2 1G  wavefunctions  are  given  in 

Table  6.  The  only  states  which  approach  the  ideal  0.75  value  (S(S+1)) 

2  2 
are  the  planar  B^  state  and  the  staggered  B^  state.  In 

these  states  the  extra  alpha  spin  electron  is  almost  completely 


localized  in  a  p  orbital  perpendicular  to  the  CHj  plane.  The  other 
states  indicate  a  large  amount  of  spin  polarization  where  alpha  and 
beta  spins  localize  in  different  regions  of  the  molecule.  This  gives 
some  indication  that  a  restricted  program  must  use  additional 
configurations  to  achieve  the  same  result. 

It  was  found  that  the  3-2 1G  optimized  geometries  are  unstable  with 
respect  to  distortions  from  The  results  are  also  presented  in 

Table  4,  5,  and  6.  The  distortions  and  energy  changes  are  larger  than 
found  for  MNDO.  When  correlation  is  used  the  distortion  becomes 
unfavorable. 

A  perturbational  correlation  treatment  such  as  the  UHP3  method 
becomes  ineffective  when  the  single  Hartree-Fock  configuration  is  not 
the  dominant  configuration  in  the  wavefunction.  Since  it  is  known 
that  there  are  several  conf igurations  of  the  same  symmetry  and  of 
comparable  energy,  a  MCSCF  treatment  was  used  to  continue  the  study. 
The  program  ALIS^  developed  by  Klaus  Ruedenberg  and  coworkers  and 
sponsored  by  the  U.S.  Department  of  Energy  was  used  to  determined 
MCSCF  energies  and  wave functions. 

In  order  to  'feel  out*  the  potential  surface  and  determine  which 

configurations  were  important  for  the  different  states  MCSCF 

18 

calculations  were  made  using  a  ST0-3G  basis  •  In  all  calculations 

IQ 

reported  below  the  3-2 1G  geometries  were  used  •  It  was  later 

confirmed  that  the  largest  contribution  to  the  wavefunction  was  the 

configuration  used  to  generate  the  geometry.  It  can  be  argued  that 

using  geometries  optimized  with  a  single  configuration  might  'favor1 

that  configuration  in  a  MCSCF  scheme,  however  the  extra  expense  of 

optimizing  at  this  level  was  not  feasible.  MCSCF  calculations  were 

2 

performed  on  the  two  staggered  A^  states  obtained  at  the  3-2 1G 

level  (vide  supra)  using  different  configurations.  For  the  geometry 
obtained  with  3  open  shells,  the  leading  contribution  (one  space 
orbital  with  two  spin  functions)  also  contained  the  3  open  shells 

(a^,  b^,  and  accounted  for  63%  of  the  wavefunction.  For  the 

geometry  obtained  with  one  open  shell  the  same  configuration 

accounted  for  49%  of  the  wavefunction. 

20 

The  steps  taken  to  obtain  the  MCSCF  energy  are  given  below. 
An  initial  guess  of  orbitals  is  made  by  doing  a  restricted  open  shell 
calculation.  For  the  ST0-3G  basis  the  energy  normally  converged  but 
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for  Che  6-31G  basis  the  energy  diverged.  In  the  latter  case  it  was 
necessary  to  stop  the  SCF  at  the  point  of  divergence  and  orbital 
improvement  was  made  at  the  expense  of  extra  MCSCF  microiterations. 
Next  a  limited  MCSCF  was  performed  with  a  few  intuitive 
configurations.  Using  this  wavefunction  a  Cl  in  a  larger  space 
(approximately  500  Spin  Adapted  Antisymmeterized  Products)  was  made  to 
determine  which  configurations  made  substantial  contributions  to  the 
wavefunction.  These  configurations  were  then  used  in  another  MCSCF 
calculation  with  the  Cl  orbitals  as  the  initial  guess.  The  process 
was  repeated  until  the  MCSCF  and  Cl  had  the  same  dominant 

configurations  and  until  no  configuration  in  the  Cl  which  was  not 
included  in  the  MCSCF  contributed  it  or  more  to  the  wavefunction.  The 
final  MCSCF  was  carried  out  with  between  5  and  11  spin  adapted 

configurations.  This  provides  an  approximation  to  the  full  500 
configuration  MCSCF  as  the  most  important  configurations  have  been 
iteratively  selected. 

The  MCSCF/ST0-3G  results  are  presented  in  Table  7.  The  first 

column  represents  the  restricted  open  shell  energies  while  the  second 

column  gives  the  MCSCF  energy  with  the  number  of  configurations  in 

parenthesis  and  the  per  cent  contribution  of  the  most  important 

configuration.  The  third  column  gives  the  Cl  energy  with  the  number 

of  spin  adapted  configurations  in  parenthesis  (8  orbitals  and  7 

electrons  in  C^v  sytmnetry).  It  has  been  pointed  out  that  in  some 

cases  carrying  out  a  calculation  in  C2v  symmetry  may  lead  to  a 

higher  energy  than  if  only  Cg  symmetry  were  used,  but  in  the 

calculations*^3^  for  N0^  it  was  found  that  the  C^v  solution  was 

the  same  as  the  Cg  solution.  For  this  reason  plus  economics,  the 

^2v  P°*nt  £rouP  was  used  when  possible. 

Relative  energies  are  presented  in  Table  9.  In  contrast  to  MNDO 

2 

and  ab  initio  results,  the  planar  A-  state  is  predicted  to  be  the 

L  2 

most  stable  of  the  C^y  states  followed  by  the  planar  B^  state 

8.0  kcal/mol  higher.  The  rotational  barrier  for  the  2B.  to 

2  1 
B2  conformational  change  is  18.6  kcal/mol,  considerablely  higher 

than  the  MNDO  (0.4  kcal/mol)  or  the  UMP3/6-31G  (6.7  kcal/mol) 

2 

barrier.  The  barrier  for  the  rotation  of  the  A^  state  (45.1 
kcal/mol)  is  larger  than  the  UMP3/6-31G  value  (51.1  kcal/mol). 


Table  7  Tot* I  Energies  for  Different  States  of 
CHjNO^  Based  on  Mul ti configurational  ST0-3G  Wavefunction 


State 

RSCF 

MCSCF 

#conf . 

%a 

Cl  #conf. 

%b 

Planar 

% 

239.7296 

239.79270 

(4) 

239.86545 

(582) 

n% 

?  39.88456 

(6) 

65% 

239.90642 

(582) 

< 

i% 

s 

239.7699 

239.86848 

(6) 

81% 

239.89358 

(588) 

< 

i% 

239.7512 

239.82680 

(5) 

87% 

239.83138 

(594) 

6% 

239.84988 

(9) 

239.85617 

(594) 

< 

1% 

Stagg. 

239.5864 

239.79005 

(8) 

52% 

239.79205 

(584) 

< 

1% 

\ 

No  SCF 

239.73403 

(7) 

239.79018 

(588) 

6% 

239.81702 

(11) 

30% 

239.82255 

(588) 

< 

1% 

\ 

239.7633 

239.85302 

(5) 

83% 

239.86399 

(594) 

< 

1% 

a)  Largest 

contribution  from  a  s 

ingle 

configuration  in 

the  MCSCF 

wave function# 

b)  Largest  contribution  from  a  single  configuration  in  the  Cl  which  was 
not  optimized  in  the  MCSCF  wavefunction* 

The  ST0-3G  minimal  basis  set  is  known  to  be  inferior  to  the  near 

2 1 

double  zeta  quality  6-31G  basis  set  •  For  this  reason  the 

calculations  were  repeated  with  the  larger  basis  set*  It  was  felt  that 

most  of  the  important  configurat ions  would  already  be  identified  from 

the  ST0-3G  basis  allowing  rapid  convergence  to  the  most  important 

configurations  in  the  expanded  basis. 

The  absolute  energies  for  the  MCSCF/6-31G  calculations  are  presented 

in  Table  8  and  the  relative  energies  are  presented  in  Table  9.  The 

restricted  open  shell  energies  are  not  given  as  the  SCF  diverged*  The 

percent  contribution  of  the  single  largest  contributing  configuration 

clearly  indicates  that  a  single  determinant  cannot  be  adequate*  The 
2  2 

planar  Bj  and  staggered  states  have  one  configuration 

comprising  approximately 

Table  8  Total  Energies  for  Different  States  of 
CHjNOj  Based  on  Mult iconf igurational  6-31G  Wavefunction 

State  MCSCF  #conf*  %'  Cl  #conf*  l ^ 

P l anar  2A?  242.86061  (6)  242.89866  (582)  less  1* 

242.9238  (8)  73%  242.94688  (582)  less  1% 


*Bi 

242.8978 

(7) 

242.90731 

(588) 

<  1% 

242.94927 

(7) 

89% 

242.95541 

(588) 

<  1% 

\ 

242.86 

(5) 

242.90560 

(594) 

<  1% 

242.92690 

(8) 

91% 

242.94086 

(594) 

<  1% 

2®1 

242.92308 

(7) 

81% 

242.92988 

(588) 

<  1% 

2ac 

2 

242.90617 

(8) 

70% 

242.92207 

(582) 

<  1% 

2A" 

242.92 

(7) 

242.92671 

(  68) 

8% 

242.93088 

(8) 

242.95703 

(1176) 

<  1% 

242.96236 

(8) 

90% 

242.97867 

(1176) 

<  1% 

Stagg.  2Aj 

242.6 

(6) 

242.62593 

(584) 

34% 

242.80496 

(7) 

242.81163 

(584) 

3% 

242.80498 

(9) 

24% 

242.81155 

(584) 

<  1% 

s 

242.8238 

(9) 

242.83152 

(584) 

<  1% 

242.82383 

(9) 

52% 

242.83165 

(584) 

<  1% 

2bi 

242.84862 

(11) 

242.86150 

(588) 

<  1% 

242.84806 

(11) 

49% 

242.86847 

(588) 

<  1% 

2b2 

242.93351 

(5) 

90% 

242.96780 

(594) 

<  1% 

a)  Largest 

contribution  from  a 

single  configuration  in  the  MCSCF 

wave function* 


b)  Largest  contribution  from  a  single  configuration  in  the  Cl  which 

was  not  optimized  in  the  HCSCF  wavefunction. 

2  2 
c)  The  B.  calculation  was  at  the  A~  geometry  and  the 

2  1  2  * 
calcuation  was  at  the  geometry (see  text). 


Table  9  Relative  Energies  (kcal/mol)  for  Different  States  of 
CH2N02  Based  on  MCSCF/ST0-3G  and  MCSCF/6-31G 


State 

MCSCF/STO-3G 

+CI 

MCSCF/6-31G 

+CI 

Planar  ^A_ 

-19.8 

-26.6 

6.1 

13.1 

s 

-9.7 

-18.6 

-9.9 

7.8 

% 

V- 

2.0 

4.9 

4.1 

-18.1 

16.9 

-6.8 

Stagg.  2Aj 

39.5 

45.1 

80.6 

98.0 

\ 

68.8 

85.4 

\ 

22.6 

46.3 

53.6 

62.3 

\ 

0 

0 

0 

0 
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90%  of  the  wavefunction.  These  corresponds  to  one  electron  in  a  p 

orbital  perpendicular  to  the  CH?  plane  and  are  the  same  states  that 
2  Z 

produced  'normal1  S  values.  Even  for  these  two  states  the  orbitals 
occupations  showed  a  significant  population  of  a  'virtual'  orbital* 
In  Figure  8  and  9  the  nonbonding  oxygen  b^  orbital  contained  1.8 

electrons  while  the  antibonding  (bonding  across  oxygen)  contained  0*2 
electrons. 

N 

This  suggests  that  there  may  be  some  stabilization  through  pi  overlap 
of  the  oxygen  p  orbitals. 

The  MCSCF/ 6-3 1G  plus  Cl  lead  to  a  predicted  staggered 

C2v  ground  state  for  the  nitromethyl  radical.  This  is  different 

from  each  of  the  previous  methods.  The  planar  ^B.  state  is  7.8 

2  1 

kcal/mol  less  stable  and  the  planar  A*  state  is  13.1  kcal/mol 

2  Z  2 
above  the  staggered  state.  The  difference  between  the 

states  is  72.3  kcal/mol  which  is  almost  identical  to  the  STO-3G 

difference  and  greater  than  that  predicted  by  UMP3/6-31G. 

The  relative  energies  have  not  converged  with  respect  to 

increasing  the  sophistication  of  the  basis  set.  Therefore  it  is 

premature  to  give  the  definitive  electronic  state  ordering.  In 

particular  polarization  functions  as  well  as  a  more  complete  Cl  may 

make  significant  contributions  to  the  relative  energies. 

Another  possible  influence  on  relative  energies  could  be 

distortions  from  C^v*  This  was  found  to  be  favorable  in  the  MNDO 

approximation  and  at  the  UHF/3-21G  and  UHF/6-31G  levels  (though  not  at 

the  UMP2/6-31G  or  UMP3/6-31G  levels).  In  order  to  see  if  the  planar 
2  2 

A^  and  states  might  interact  upon  loss  of  the  sigma  (xz) 

plane  (vide  supra) ,  MCSCF  calculations  were  carried  out  on  the  planar 

2  2  2 
A^  state  at  the  B^  geometry  and  on  the  B^  state  at  the 

2A^  geometry.  The  results  shown  in  Figure  10  (see  Table  8  for 

absolute  energies)  indicate  that  a  crossing  must  occur  somewhere 

between  the  two  geometries  along  a  reaction  coordinate* 

A  MCSCF  calculation  was  carried  out  on  the  distorted  C2v 


geometry  optimized  with  the  UHF/3-21G  basis  (see  Table  6).  The 

results  given  in  Table  8  and  9  indicate  that  the  distortion  is 

2 

favorable  and  yields  a  A"  state  14.6  kcal/mol  more  stable  than  the 
2 

planar  B.  state*  This  is  6*8  kcal/mol  more  stable  than  the 

l2 

staggered  and  thus  becomes  the  ground  state*  It  would  be 

important  to  precisely  determine  the  Cg  minimum  before  structural 
predictions  are  made  as  the  UHF/3-21G  optimized  geometry  may  not  be 
accurate  %  but  it  is  felt  that  such  optimization  would  not 

qualitatively  change  this  result* 


calculations  at  the  four  labeled  points* 
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RECOMMENDATIONS 


On  the  basis  of  the  above  study  the  following  recommendations  can 
be  made: 

(1)  That  single  determinant  methods  for  the  nitromethyl  radical  may 
neglect  important  contributions  to  the  wavefunction  and  must  be  used 
with  some  caution  in  making  structural  and  energetic  predictions. 

(2)  That  semi-empirical  and  ab  initio  methods  provide  similar 
qualitative  predictions,  however  it  is  more  difficult  to  obtain 
different  electronic  states  with  the  MNDO  method. 

(3)  That  distortions  from  symmetry  must  be  considered  when 

states  differing  by  only  one  spin  orbital  cross  at  some  geometry. 
(A)  That  unrestricted  ab  initio  results  without  correlation  probably 
overestimate  such  distortions. 

(5)  That  a  similar  study  be  undertaken  at  the  ab  initio  level  for  the 
radicals  derived  from  1-nitropropene  which  may  be  a  more  realistic 
model  for  the  radicals  produced  from  TNT. 

(6)  That  in  order  to  study  the  decomposition  mechanism  for 
nitromethane  the  reaction  pathway  from  the  nitromethyl  radical  to 
CH2  plus  N(>2  end  to  CI^O  plus  NO  be  studied. 
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